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Abstract 

Theories of interacting gauge fields and fermions can possess a running gauge couphng with an 
infrared attractive fixed point (IRFP). We present a minimal description of the physics of these 
systems and comment on some simple expectations for results from lattice simulations done within 
the basin of attraction of the IRFP in these theories. 
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I. INTRODUCTION 



Interest in non-Abelian lattice gauge theories which differ from QCD in having more 
quark flavors or quarks in larger representations has been growing in recent years. These 
studies aim to address physics that may appear at the energy of the Large Hadron Collider 
or beyond, described by theoretical constructs such as technicolor [l[ or "unparticles" [2|. 

Most of the interesting theories are asymptotically free and an infinite cut-off limit can 
be defined at the ultraviolet Gaussian fixed point (FP) at bare couplings g'^ = 0, m = 
0. Starting with bare parameters near the FP both the running gauge coupling and the 
mass increase as the energy scale is lowered. As long as the coupling g"^ remains small 
the perturbative renormalization group f3 function describes its running. There are two 
possibilities at stronger coupling [sj, 0]. The first is when the Gaussian fixed point is the 
only fixed point with a diverging correlation length. A familiar example is QCD with two 
light flavors. The f3 function stays negative even at strong coupling, chiral symmetry is 
spontaneously broken and the model is confining. The second possibility is that the (3 
function develops another zero, corresponding to an infrared attractive fixed point (IRFP). 
In this case the infrared physics of the massless theory will possess a conformal symmetry 
that precludes confinement and the spontaneous breaking of chiral symmetry. These are 
the "unparticle" theories. Just before the conformal window opens up there is a possibility 
that the j3 function, while staying negative, becomes small, thus the running of the coupling 
slows down to "walking". Such theories are leading candidates for technicolor models. 

Most of the interesting cases in the literature have nonperturbative dynamics. The the- 
oretical approaches used to analyze their behavior are not under complete control. Accord- 
ingly, several groups have begun to use lattice methods to investigate these models. In 
at least three cases (an SU{3) gauge group with 12 flavors of fundamental representation 
fermionsjl, [gI, SU{3) with two flavors of symmetric-representation fermionsQ], and SU{2) 
with two flavors of adjoint representation fermionsjsl) authors have reported evidence for an 
IRFP theory. 

This short note is written to present the simplest theoretical description of IRFP theories. 
While we believe that a considerable fraction of the lattice community is familiar with the 
physics we will describe, we have participated in enough conversations and read enough 
published papers in which the description we will give was incompletely presented, that 
we feel a review might be useful. A large part of the non-lattice literature about IRFP 
theories is devoted to describing the transition from an IRFP system to a confining one, as 
the number of flavors of fermions is reduced [ol, [l^, 11, l^, 13|. Our "theoretical minimum" 



might be useful as a benchmark against which simulation results and these predictions can 
be compared. 

After some introductory remarks, we would like to make two points: 

• When a theory is in the basin of attraction of an IRFP, its gauge coupling is irrelevant. 
The only relevant coupling, which controls the leading scaling behavior of the theory, 
is the fermion mass. 

• The value of the IRFP gauge coupling, g*, is scheme- dependent. Spectral observables 
cannot depend on its value, and can only be sensitive to the extreme values of the 
bare gauge coupling which mark the basin of attraction of the IRFP region. The 
renormalization group flow is always towards g*, but its absolute direction depends on 
the scheme- dependent location of g*. 
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We now elaborate on these points and their consequences. 



II. MAPPING THE PHASE STRUCTURE WITH THE RENORMALIZATION 
GROUP 

The renormahzation group (RG) transformation is a frequently used tool in studying the 
infinite cut-off limit of quantum field theories. There are two different approaches. The 
Callan-Symanzik equation describes the problem from the perturbative point of view and 
is most often used in connection with continuum regularization schemes. It describes the 
change of the parameters of the theory as the function of the cut-off, or equivalently, the 
renormahzation subtraction mass /i. The calculation is perturbative and only parameters 
already present in the Lagrangian are considered. For example the /3 function 

describes the running of the coupling gf^, and there are similar expressions for the mass 
and any other couplings present. Zeros of the (3 function correspond to either ultraviolet 
or infrared fixed points (UVFP and IRFP), depending on the slope of /3{g'^) (negative or 
positive, respectively). 

The other option is the inherently non-perturbative Wilson RG approach. There one con- 
siders the evolution of all the possible couplings of the system under an RG transformation 
that preserves the internal symmetries of the system but integrates out the cut-off level UV 
modes. The fixed points of the transformation are characterized by the number of relevant 
couplings or operators, i.e. couplings that flow away from the FP. Continuum (or infinite 
cut-off) limits are defined when the relevant couplings are tuned towards the FP. Irrelevant 
couphngs flow in the IR to values which are independent of their UV values. The number of 
relevant operators and their speed along the RG flow lines arc universal properties, related 
to the critical properties of the underlying continuum limit. However the location of the FP 
is not physical, in fact different RG transformations have different fixed points. 

The relation between the Callan-Symanzik and Wilson RG approaches are straightfor- 
ward: the Callan-Symanzik P function describes the RG flow of one or a few couphngs, 
usually implicitly in perturbation theory. The zeroes of the f3 function correspond to the 
fixed points of the Wilson RG. While the Callan-Symanzik equations can provide the first 
hint about the phase structure, unless the predictions are well controlled by perturbation 
theory, a non-perturbative approach is needed to verify the existence and study the proper- 
ties of the relevant fixed points. Lattice calculations provide an excellent tool for that. We 
believe that the Wilson RG approach is better suited to interpret lattice results and in the 
following we will use mainly that language. 



III. THE PHASE DIAGRAM OF GAUGE THEORIES WITH MANY FERMIONS 

Consider a theory with SU (Nc) gauge group, coupled to Nf flavors of fermions in rep- 
resentation R. In our mind is a lattice simulation and that is the phase diagram we will 
describe. The Gaussian fixed point at = 0, m = is well understood perturbatively. 
The mass is a relevant operator and it will presumably will remain so even at strong gauge 
coupling, m = is a critical surface and we want to investigate the running of the gauge 
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coupling along it. In the perturbative region this is described by the (3 function which is 
universal up to two loop level: 



(3) 



If the fermion number A'^^ < {ll/A)Nc/T{R) (16.5 for fundamental representation fermions 
in QCD), 61 < and both the gauge coupling and the mass are relevant operators at the 
Gaussian FP. The theory is asymptotically free, the coupling g"^ increases as the lattice 
spacing (or inverse cut-off) increases. For small number of fermions the perturbative (3 
function remains negative even at strong coupling. Lattice simulations with Nf = 2, 2+1, 
and to some extent up to Nf < 8 verify this: QCD with 2-8 light fermions is confining and 
chirally broken everywhere . With increasing Nf the two-loop term in Eq. [2] changes 
sign, suggesting a possible zero in f3{g). When the quarks are in higher representations of 
the gauge group, this happens more easily. The zero of l3{g) corresponds to an IRFP in the 
gauge coupling. That is, the gauge coupling is irrelevant, and the only relevant parameter 
is the mass term. The continuum limit is defined in the basin of attraction of this IRFP by 
tuning the mass to zero. 

The existence of the perturbatively predicted IRFP and the properties of the theory at 
even stronger coupling have to be studied with non-perturbative methods, like lattice simu- 
lations. The lattice provides an UV regularization scheme which is equivalent to continuum 
schemes at weak coupling. At strong gauge coupling lattice artifacts can completely change 
the system. It is generally believed that in the strong cou plin g limit lattice models are 



always confining and chirally broken, independent of Nf |l5l. Il6l. Il7l|, This was observed in 



early Nf = 16 simulations |18l. Il9l|. though numerical results in Ref. [20| contradict this. In 
any case, the details in the strong coupling region are not universal and, for now, not very 
important either. Fig. [1] shows a conjectured phase diagram in the g^ vs m plane when an 
IRFP exists. The thick line at large g"^ shows the (possible) confining and chirally broken 
phase, but around the IRFP the theory is chirally symmetric and deconfined. The dashed 
line indicates the phase boundary between the two phases. It is likely only a crossover at 
finite mass. We note that in a recent paper it was suggested that in these many-fermion 
theories with an IRFP there is always another FP (zero of the (3 function) at a stronger 



coupling [13|. If this conjecture is valid, the phase diagram has another UVFP at m = 0. 
This could be the boundary between the conformal and strong coupling phases in Fig. [H or 
an entirely new critical point. 

It is important to note that while at a UVFP an observable phase transition occurs, there 
is nothing observable at the IRFP. The location of an IRFP on the critical surface is not 
physical. It depends on the specific renormalization group transformation, and therefore no 
physical observable can identify it. In particular, this means that studies of candidate theo- 
ries which only measure spectroscopic observables cannot directly detect coupling constant 
flow through the dependence of observables on irrelevant bare couplings. They can only 
detect the perimeter of the basin of attraction of the FP. 
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FIG. 1: A conjectured phase diagram for a lattice theory with an IRFP. The thick hne at m = is 
the confining, strong-coupling region of the massless theory. The dotted line separates a confining 
phase (which includes the entire line at large fermion mass) from the basin of attraction of the 
IRFP. Arrows show some IR flow lines and the renormalized trajectory extending from the critical 
point. 

IV. ANALYSIS OF THE PHASE DIAGRAM 

To begin, note first that the gauge coupling plays a very different role around the UVFP 
and in the vicinity of the IRFP or trivial phases: at the UVFP g is relevant, and to drive 
the correlation length(s) to infinity in units of the cutoff scale involves a fine tuning of g 
toward zero. (To be precise, the gauge coupling is marginally relevant with respect to the 
Gaussian fixed point, because Pi^g"^) ~ big^.) In the other phases, IR physics is independent 
of the value of g at the UV scale, because the IR flow of the gauge couphng is into its 
flxed point value: the gauge coupling is irrelevant. This is a simple example illustrating 
the general feature, that relevance or irrelevance of an operator is measured with respect 
to a fixed point. The presence of the IRFP means that the linearized behavior of the beta 
function near the fixed point is 



where s is the (IR) scale change, /i e~*/i, and ^/^ < is the scaling exponent of the gauge 
coupling. The running coupling then is g'^{s) = g^* + {gi — g^*)^-'^^^ , i-e. g{s) — >• 51* as s — >• oo. 
We discuss the value of yg below. 

Consider next the general situation, a theory in D Euclidean dimensions defined by a 
set of scaling operators {u}. Under a real-space blocking by a factor b, the operators run 
multiplicatively to/from their FP values, u[ — Uio = W^{ui — Uio), with the usual assignment 
of relevancy, irrelevancy, or marginality, depending on whether yi is positive, negative, or 
zero. Let us assume that the mass is the only relevant coupling and its FP value is zero; its 
exponent we will label as y^. Only the leading exponent governs the correlation length ^, 




(4) 
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which diverges with the usual algebraic behavior 



^ ~ m . (5) 

A standard textbook analysis tells us that the the singular part of the free energy per site (it 
is the part of the free energy containing all the non-analyticity at the critical point) scales 
as 

/s(m, Ui) = m^'^"^ fsiniQ, Um + {ui - Uio) I — ), (6) 

which can be Taylor expanded as 

fs{m) = m^/^'"(Ai + Asm'^'l/s^'-), (7) 

where Ai and A2 are non-universal constants. All other observables have a similar expansion. 

Precisely at m = the entire basin of attraction of the FP is critical - at long distance 
all correlation functions decay algebraically 

{Mr)<Pm) = j:^^. (8) 
j 

This behavior is only achieved by setting the quark mass to zero. The nonanalytic depen- 
dence of the free energy on the fermion mass is responsible for the algebraic behavior of 
correlation functions in the massless theory. Eq. [8] is only true asymptotically. At distances 
which are comparable to the UV cutoff (lattice spacing a) there are additional contributions. 
They arise from the non-singular part of the free energy, from physics at the intermediate 
scales which has been integrated out in the construction of Eq. O These give extra contri- 
butions going as exp(— r/a). 

On the m = critical surface all couplings are irrelevant; they all flow into the IRFP. 
The interesting physical quantity near the FP is the leading relevant exponent. This is not 
accessible from a lattice-based RG study at m = (for example, a conventional Schrodinger 
functional study which computes P{g)). These studies only give yg because by construction 
they are only sensitive to flow into the FP. They tell us that there is an IRFP, no more. 

Unitarity bounds for D = A conformal field theories [2l| constrain the scaling dimension 
of the condensate (ipip) ~ /i"^ to lie in the range 3 > 7 > 1. Since mipilj is scale invariant 
and since all dimensionful quantities scale with the correlation length as in Eq. 7 = 
4 — 1/m- (There are other definitions for the scaling dimension in the literature but this one 
is consistent with the bound of Ref. [21].) We do not want to discuss specific techniques 
for determining it in a simulation, and defer a discussion of possible values for i/m to later 



work j22|. In the meantime, however, there is one case for which it is easy to expose the 
operator scaling hierarchy [i, 23]. This is a lattice theory of Nc colors and Nf flavors in some 
representation, and that the limit Nc 00 and Nf/Nc fixed, so that the Banks-Zaks FP 
occurs at a small value of g* ~ e. Engineering dimensions give basically everything. For this 
theory, the condensate has scaling exponent 7 = 3 — 0(e). The gauge coupling contributes a 
scahng exponent l/c, ~ e (which can be computed in perturbation theory), and lattice-based 
artifacts contribute negative ?/j's. 

This is the upper end of the conformal window. What happens at the lower end of the 
conformal window, where the system converts from being an IRFP theory to a confining, 
chirally broken theory, is a subject of long-standing interest. Observations of ym will likely 
be needed to elucidate this point. 
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V. CONSEQUENCES FOR SIMULATIONS 



One feature of IRFP theories is concealed by the Callan-Symanzik formahsm: the location 
of a fixed point is not universal. Different choices of renormalization group transformations 
will result in different values of g* . This has several consequences for simulations. 

First of all, observables taken from correlation functions cannot depend on the actual 
value of the FP coupling determined in some scheme: the correlation function does not 
know about the scheme. The important part of expressions like Eq. [7] is the statement that 



Real QCD (for example, an S'[/(3) gauge theory with a small number of flavors of fun- 
damental representation fermions) is more complicated than an IRFP theory. Now there 
are two relevant couplings which flow away from the {g'^,mq) = (0,0) fixed point. Lattice 
observables (like masses) depend strongly on both of these couplings. Indeed, making lattice 
predictions typically involves first finding a map of mq{g'^) along which some mass ratio is 
fixed, then taking the continuum limit by moving along this line toward the UVFP. If one 
tried to do this for an IRFP theory, one would be trying to express (or tune) the relevant 
coupling in terms of an irrelevant coupling. 

Next, because different choices for renormalization schemes result in different values for 
g*, it is easy to imagine a situation where the running coupling in one scheme is increasing 
under flow to the IR, while the running coupling in another scheme is decreasing. In fact, 
this is another diagnostic for the presence of an IRFP: one can vary the renormalization 
group transformation and see if the direction of coupling constant flow changes. 

Much continuum beyond-Standard Model phenomenology uses a scale dependent coupling 
constant defined in some particular scheme, such as MS. A lattice calculation of such a 
coupling involves both a lattice simulation and an additional scheme matching calculation. 
Critital exponents such scheme-independent . 

Simulations are carried out at particular values of the bare couplings. Eq. [7] tells us that 
the influence of the gauge coupling on observables dies away in the zero mass limit. In a 
sense, varying the gauge coupling while remaining near the basin of attraction of the IRFP 
is like doing lattice simulations for ordinary QCD with several choices of the lattice action 
- in the scaling limit, the answer is not supposed to depend on the choice of action. Away 
from infinite correlation length, different lattice actions will give different predictions for 
observables. But these differences are just scale violations. In this case, the choice of bare 
gauge coupling is as much a choice of the action as a particular discretization would be. 

The correlation length ^ only varies with the bare mass according to Eq. Owhen the system 
size L is much larger than S,- This is probably the the most serious practical constraint on 
the choice of bare parameters. 

The size of the exponent Ug also acts to minimize the dependence of simulation results 
on g^: It is likely the value of Ug is very small. The exponent can be seen in simulations. 
The authors of Ref. [1] report Ug ~ —0.15. The authors of Ref. [1] do not quote a number, 
but they give a plot of the j3 function from which a value of about -0.16 can be inferred. 
The authors of Ref. Q only show an integrated beta function, but parameterizing their data 
with a simple model (a second order formula with the perturbative hi and a 62 chosen to 
give a zero) also gives a small yg ~ —0.06. Gardi and Grunberg(l2| present perturbative 
calculations for both QCD-like and supersymmetric theories. Only at the bottom of the 
conformal window, and only for non-supersymmetric theories, is —yg greater than about 
0.5. If that is the case, then in formulas like Eq. [7] the correction term is nearly a constant. 



the size of scale violations in observables is 
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This has two consequences: 

• In a lattice simulation of a particular system which is a candidate for being an IRFP 
theory, the observation that long distance observables become independent of g at 
small m is a positive indication of IRFP behavior 

• In the basin of attraction of an IRFP, one does not need to fine tune the gauge coupling; 
one can compute wherever it is convenient. 

What evidence is there for the behavior we have described, from published results for 
candidate IRFP theories? Two groups have studied the SU{2) gauge group with two flavors 
of adjoint fermions. This model has a confining strong coupling phase and a weak coup ling 
phase which the authors of Ref. [1] have argued is conformal. The authors of Ref. [2^ 
measure the pseudoscalar and vector masses in finite volume at zero quark mass. Finite 
volume limits the correlation length to be proportional to the system size. They observe 
(in their Fig. 5) that these masses are independent of the bare gauge coupling in the weak 
coupling phase. They also note the constancy of Creutz ratios with respect to gauge coupling 



(ratios of Wilson loops) in this phase; see their Fig. 7. The authors of Ref. [25| observe 
that their pseudoscalar and vector masses are independent of gauge coupling in the weak 
coupling phase; compare their figures 5-8. 

One of us collaborated with Svetitsky and Shamir in studies of SU (3) gauge theory with 
two flavors of sextet fermions 2^. This system also has a weak coupling phase which is a 



candidate for an IRFP phase. It is easier to replot data published there: Fig. [2] shows the 
pseudoscalar mass as a function of the quark mass, from simulations performed at many 
values of the gauge coupling (/? = Q/ g^) in the weak coupling phase. The finite simulation 
volume prevents the vanishing of the pseudoscalar mass at zero quark mass. Results are 
similar to what is seen in the SU{2), adjoint case, i.e., weak dependence of the meson mass 
on the gauge coupling. 



VI. CONCLUSIONS 

We have presented a minimal scenario for physics within the basin of attraction of an 
IRFP. There is one relevant coupling, the fermion mass, which governs the correlation length 
through Eq. [5l All other couplings are irrelevant. This scenario gives a context for inter- 
preting simulation results and has observable consequences. Probably, convincing evidence 
for an IRFP requires RG studies. However, the simple behavior we have described provides 
additional markers for the description of an IRFP theory with a single relevant operator. 
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